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Abstract 

Ph 

I We study the determination of the symmetry that stabilizes a dark matter (DM) candidate produced at 

P-i 

colliders. Our question is motivated per se, and by several alternative symmetries that appear in models 
r^. that provide a DM particle. To this end, we devise a strategy to determine whether a heavy mother particle 

decays into one visible massless particle and one or two DM particles. The counting of DM particles in 
these decays is relevant to distinguish the minimal choice of Z2 from a Z3 stabilization symmetry, under 
which the heavy particle and the DM are charged and the visible particle is not. Our method is novel in 
that it chiefly uses the peak of the energy spectrum of the visible particle and only secondarily uses the 
Mt2 endpoint of events in which the heavy mother particles are pair-produced. We present new theoretical 
results concerning the energy distribution of the decay products of a three-body decay, which are crucial for 
^""^ our method. To demonstrate the feasibility of our method in investigating the stabilization symmetry, we 

. 1^ apply it in distinguishing the decay of a bottom quark partner into a b quark and one or two DM particles. 

The method can be applied generally to distinguish two- and three-body decays, irrespective of DM. 
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1 Introduction 

Extensions to the Standard Model (SM) of particle physics are motivated for various reasons; per- 
haps the most important among these is the necessity of a fundamental mechanism for electroweak 
symmetry breaking (EWSB). Additionally, the related Planck- weak hierarchy problem of the SM 
must also be addressed. In such extensions of the SM, there generally exists a new particle at or 
below the TeV scale which cancels the quadratic divergence of the Higgs mass from the top quark 
loop in the SM. Such a particle is typically a color triplet with a significant coupling to the SM top 
quark, and has an electric charge of +2/3. Following the literature, we will generically call such 
particles "top partners" and denote them by T' ^ . These top partners often come along with bottom 
partners, which we similarly denote as B' . The typical reason for this is that the left-handed (LH) 
top quark is in a doublet of SU (2) l with the LH bottom quark. We then expect top and bottom 
quark-rich events from the production and decay of these new particles at the LHC. 

Another seemingly unrelated motivation for new physics at the TeV scale is the evidence for 
the existence of dark matter (DM) in the Universe, combined with the absence of a viable DM 
candidate in the SM [1]. A well- motivated candidate for this DM is found in a stable weakly 
interacting massive particle (WIMP), especially one that arises as part of an extension to the SM 
at the TeV scale. The motivation for this new physics becomes even stronger when the extension 
to the SM solves other problems inherent in the SM. These scenarios often involve heavier new 
particles that are charged under both the symmetry that keeps the DM stable and the SM gauge 
group. These new particles should then be copiously produced at the LHC and must decay into 
DM particles and SM states, given that the latter are not charged under the DM stabilization 
symmetry. Thus we expect this new physics to give rise to events at the LHC with large missing 

^In our work this name applies as long as the partners have interactions with the relevant SM particle, even if the 
partners do not directly cancel the Higgs mass divergence. 
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energy, in association with jets, leptons, and photons. 

Combining the above two hnes of argument, we reahze that the most attractive scenarios are 
those extensions of the SM which not only solve the Planck-weak hierarchy problem, but also have 
a WIMP DM candidate. In this case, it is likely that the top and bottom partners are also charged 
under the DM stabilization symmetry. These extensions will then result in top and bottom quark- 
rich events at the LHC in which the new particles give rise to missing energy. The classic example 
of such an extension is SUSY, where i?-parity stabilizes the DM [2]. The associated signals from 
the scalar top and bottom partners have been studied in great detail. A more recent example is 
little Higgs models [3] with T-parity [4]. Like SUSY, the signals from the fermionic partners of the 
top and other quarks in these models have been thoroughly studied. In short, we find that a search 
for events with top or bottom quarks and missing energy should be a top priority of the LHC. 

Once the existence of new physics has been established, the most urgent issue that will then have 
to be addressed is the determination of the details of the dynamics underlying this new physics. 
In particular, it will be crucial to determine the properties of the top and bottom partners using 
as model-independent an approach as possible. This detailed study would also offer major hints 
regarding the resolution of the Planck-weak hierarchy problem. For largely model-independent 
work on fermionic bottom and top partners' discovery potential at the LHC see Refs. [o, (>] and for 
the determination of generic partners' spin and mass see Refs. [7]. 

However, we remark that in these works it has been assumed that the top or bottom partner 
decays into only one DM particle, which is expected when the DM is stabilized by a Z2 symmetry. 
While Z2 is perhaps the simplest DM stabilization symmetry, it is by no means the only possibility: 
see references [8, 9]. The point, especially in the case of such non-Z2 symmetries, is that more than 
one DM can appear in the decays of top and bottom (and other SM) partners: for example, two 
DM are allowed with Z3 as in [S], but not with Z2. 

We believe that a truly model-independent approach to the determination of the top and bottom 
partners' properties should include this possibility of multiple DM in addition to different spins for 
the top and bottom partners. With this goal in mind, we aim to devise a strategy that uses 
experimental data to determine the number of DM in these decays and accordingly to identify the 
stabilization symmetry of the dark matter. Below, we outline a general strategy and then apply it 
to the specific case of bottom partner decays. 

We concentrate on the distinction between two general decay topologies: 

A^bX and A^bXY (1) 

where 6 is a (single) SM visible particle, X and Y are two potentially different invisible particles 
and j4 is a heavier particle that belongs to the new physics sector. In the context of the models 
that we have discussed, A is the heavy particle charged under the DM stabilization symmetry and 
the particles labeled X and Y are the DM particles. In particular, we focus on scenarios where the 
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two decays are mutually exclusive, i.e. where the stabilization symmetry and the charges of the 
involved particles are such that one decay can happen and not the other. This mutual exclusivity 
can be the case with both Z2 and Z3 as the stabilization symmetry. To wit, if the SM particle b 
is not charged under the stabilization symmetry and all the new particles A, X, Y are, then the Z2 
symmetry allows only for two-body decays of A. On the other hand, both the two and three-body 
decays of A are allowed by the Z3 symmetry by itself. However, we assume that other considerations 
forbid (or suppress) the two-body decay in this model. We choose to concentrate on this realization 
of the Z3 -symmetric model in part because this is the case that cannot be resolved using the results 
of previous work on the DM stabilization symmetry. This is the case, for instance, in Ref. [10], 
where purely two-body decays of A could be distinguished from mixed two- and three-body decays, 
but not from the purely three-body decays that we are now taking into consideration. 

In this paper, we develop a method based primarily on the features of the energy distribution 
of the visible final state b to differentiate between the cases of purely two- and three-body decays. 
We remark that this is the first work to use the energy distribution of the the decay products to 
study the stabilization symmetry of the DM. In fact, other work has typically focused on using 
Lorentz invariant quantities or quantities that are invariant under boosts along the beam direction 
of the collider. This is the case for Refs. [10, 11, 12, 13]. In particular, Refs. [10, 11, 12] used the 
endpoints of kinematic distributions to probe the stabilization symmetry of the DM, whereas our 
method relies quite directly on peak measurements and only marginally on endpoint measurements. 
Additionally, we note that the methods developed in Refs. [12, 13] apply only to the case where 
there are more than one visible particle per decay. Therefore, our result for cases where there is 
only one visible particle per decay is complementary to the results of the above references. 

Our basic strategy is explained in the following. It relies on a new result: assuming massless 
visible decay products and the unpolarized production of the mother particles, we will show that in 
a three-body decay the peak of the observed energy of a massless decay product is smaller than its 
maximum energy in the rest frame of the mother. This observation can be used in conjunction with 
a previously observed kinematic characteristic of the two-body decay to distinguish the stabilization 
symmetry of the DM. Specifically, it was shown in Ref. [14, 15] that for an unpolarized mother 
particle, the peak of the laboratory frame energy distribution of a massless daughter from a two- 
body decay coincides with its (fixed) energy in the resi-frame of the mother. 

Clearly, to make use of these observations in distinguishing two from three-body decays, we need 
to measure the "reference" values of the energy that are involved in these comparisons. Moreover, 
the procedure that is to be used to obtain this reference value from the experimental data should be 
applicable to both two and three-body decays. To this end, we find that when the mother particles 
are pair produced, as happens in hadronic collisions, the Mt2 variable can be used. Thus, these 
observations make counting the number of invisible decay products possible by looking only at the 
properties of the single detectable particle produced in the decay. However, it is worth noting that 
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our proof of the above assertion regarding the kinematics of two- and three-body decays is only 
vahd with a massless visible daughter and an unpolarized mother. Therefore, care must be taken 
when discussing cases with a massive daughter or a polarized mother. 

To illustrate the proposed technique, we will study how to distinguish between pair-produced 
bottom partners each decaying into a b quark and one DM from pair-produced bottom partners each 
decaying into a b quark and two DM particles at the LHC ^. As discussed above, a bottom partner 
appears in many motivated extensions to the SM, so we posit that this is a relevant example. 
Furthermore, we remark that the b quark is relatively light compared to the expected mass of 
the bottom partner, so that our theoretical observation for massless visible particles is expected 
to apply. Additionally, the production of bottom partners proceeds dominantly via QCD and is 
thus unpolarized. In this sense, the example of a bottom partner is well-suited to illustrate our 
technique. Finally, it is known that the backgrounds to the production of bottom partners may be 
rendered more easily manageable than for those of top partners [5], which would be a well- motivated 
alternative example. 

Specializing to the example of bottom partners, our goal then is to distinguish the two processes 
illustrated in Figure 1 at the collider 

pp B'B' bbxx for Z2 , (2) 

pp B'B' bbxxXX for , (3) 

where x is an invisible particle and a bar denote anti-particles. In these processes, we assume 
that there are no on-shell intermediate states. We consider the case where the decay into two 
X can happen only if the stabilization symmetry of the DM is Z3, while the decay into one x 
is characteristic of the Z2 case. As said before, we focus on this scenario because it has thus far 
been left uninvestigated by previous studies on the experimental determination of the stabilization 
symmetry of the dark matter [10, 12]. 

From here, we organize our findings as follows: In Section 2, we review the current theory 
and we derive new results about the energy spectrum of the decay products of two- and three- 
body decays. These are then the foundation of the general technique presented in Section 3 for 
differentiating decays into one DM particle from those into two DM particles. In Section 4, we 
apply this technique to the specific case of bottom partners at the LHC. We conclude in Section 5. 

2 Theoretical observations on kinematics 

We begin first by reviewing the relevant theoretical observations about the kinematics of two- 
body and three-body decays. Specifically, we review the remarks on two-body decays described 

^To the best of our knowledge, none of the earher work on distinguishing DM stabihzation symmetries at coUiders 
has studied this specific case. 



5 



Figure 1: The signal processes of interest for TLi (left panel) and Z3 (right panel) stabilization 
symmetry of the dark matter particle x. 



in [14]. We then generalize this result to three-body decay kinematics and study the features 
that distinguish it from two-body decay kinematics. We also briefly review applications of the 
kinematic variable Mt2 to two-body and three-body decays and discuss the distinct features of the 
two different decay processes [10, 16]. 

For the two-body decay, we assume that a heavy particle A decays into a massless visible 
daughter h and another daughter X which can be massive and invisible: 

A^hX. (4) 

On the other hand, for a three-body decay the heavy particle A decays into particles 6, X and 
another particle Y 

A^hXY . (5) 

Like particle X, particle Y can also be massive and invisible, but it is not necessarily the same 
species as particle X. 

2.1 The peak of the energy distribution of a visible daughter 
2.1.1 Tvi^o-body decay 

It is well-known that the energy of particle h in the rest frame of its mother particle A is fixed, 
which implies a b function-like distribution, and the simple analytic expression for this energy can 
be written in terms of the two mass parameters and mx ■ 

2 2 

Typically, the mother particle is produced in the laboratory frame at colliders with a boost that 
varies with each event. Since the energy is not an invariant quantity, it is clear that the 5 function- 
like distribution for the energy as described in the rest frame of the mother is smeared as we go 



'We will consider the case of a massive daughter particle in future work. 
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to the laboratory frame. Thus, naively it seems that the information encoded in eq. (6) might be 
lost or at least not easily accessed in the laboratory frame. Nevertheless, it turns out that such 
information is retained. We denote the energy of the visible particle b as measured in the laboratory 
frame as Eh. Remarkably, the location of the peak of the laboratory frame energy distribution is 
the same as the fixed rest-frame energy given in eq. (6): 

^peak^^*^ (7) 

as was shown in [14, 15]. 

Let us briefly review the proof of this result while looking ahead to the discussion of the three- 
body case. As mentioned before, the rest-frame energy of particle b must be Lorentz-transformed. 
The energy in the laboratory frame is given by 

Eb = E;j{1 + /? cos 9*) = El{-i + - 1 cos 6*) , (8) 

where 7 is the Lorentz boost factor of the mother in the laboratory frame and 0* defines the angle 
between the emission direction of the particle b in the rest frame of the mother and the direction 
of the boost /3, and where we have used the relationship 7/3 = \/'^'^ — \. If the mother particle is 
produced unpolarized, i.e., it is either a scalar particle or a particle with spin produced with equal 
likelihood in all possible polarization states, the probability distribution of cosO* is flat, and thus so 
is that of Eb. Since cos 9* varies between —1 and +1 for any given 7, the shape of the distribution 
in Eb is simply given by a rectangle spanning the range 



Eb G 



El (7 - x/^^) , El (7 + x/^^)] . (9) 



It is crucial to note that the lower and upper bounds of the above-given range are always smaller 
and greater, respectively, than Eb = El for any given 7, so that El is covered by every single 
rectangle. As long as the distribution of the mother particle boost is non-vanishing in a small 
region near 7 = 1, i?* is the only value of Eb to have this feature. Furthermore, because the energy 
distribution is flat for any boost factor 7, no other energy value has a larger contribution to the 
distribution than El. Thus, the peak in the energy distribution of particle b is unambiguously 
located at Eb = El. 

The existence of this peak can be understood formally. From the fact that the differential decay 
width in cos^* is constant, we can derive the differential decay width in Eb for a fixed 7 as follows: 



1 dT 
fdEb 



1 dr dcosO* 



fixed 7 rdcos^* dEb 



fixed 7 



1 



2Ely^i^ 



e 



Eb - El (7 - V7'-l)l e \-Eb + El (7 + Vt'-i)! (10) 



where the two Q{Eb) are the usual Heaviside step functions, which here merely define the range 
of Eb. To obtain the full expression for any given Eb, one should integrate over all 7 factors 
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contributing to this Eh. Letting ^'(7) denote the probabihty distribution of the boost factor 7 
of the mother particles, the normalized energy distribution /2-body(-£'fe) can be expressed as the 
following integral 

9(7) 



00 



A^.AE.)-J^^^^^^^^^. (11) 

The lower limit in the integral can be computed by solving the following equation for 7: 

Eh = E;(j±y'j^-l) (12) 



with the positive (negative) signature being relevant for Eb > E^ {Eb < E^). We can also calculate 
the first derivative of eq. (11) with respect to Eb as follows: 

^b \ I I , -t^b 



The solutions of f2-hody(Eb) = give the extrema of f2-hody{Eb), and given the expression 
fLhodyi^b) in eq. (13), these zeros originate from those of 5(7). For practical purposes, one can 
take 5(7) to be non- vanishing for particles produced at colliders for any finite value of 7 greater than 
1 ^. As far as zeros are concerned, two possible cases arise for g{l) (corresponding to Eb = E^). If it 
vanishes, f'2-hody^-Eb = El) oc g{l) = 0, which implies that the distribution has a unique extremum 
at Eb = El If g{l) / 0, 

body(-^'') overall sign change at Eb — E^ As a result, the 

distribution has a cusp and is concaiV6-d-Own at Ei) — Moreover, tlie function ^2-body 

(Eb) has 

to be positive to be physical, and has to vanish as Eb approaches either or oo, which is manifest 
from the fact that in those two limits the definite integral in eq. (11) is trivial. Combining all of 
these considerations, one can easily see that the point Eb = El is necessarily the peak value of the 
distribution in both cases. 

2.1.2 Three-body decay 

We now generalize the above argument to three-body decays. We denote the energy of the visible 
particle b measured in the rest frame of the mother particle A as Eb- We also denote the normalized 
rest-frame energy distribution of particle b as h{Eb). In the two-body decay, this rest-frame energy 
is single- valued (see eq. (6)), and so the corresponding distribution h{Eb) was trivially given by 
a J-function. However, when another decay product is introduced, for instance, particle Y in 
eq. (5), then the energy of particle b is no longer fixed, even in the mother's rest frame: h(Eb) ^ 
5 (^Eb — El) . Although the detailed shape of this rest-frame energy distribution is model-dependent, 
the kinematic upper and lower endpoints are model- independent. Since particle b is assumed 



''it must be noted that due to the finite energy of the coliider, there is a fcinematic upper hmit for the boost factor 
7 of the heavy mother particles. However, this kinematic limit is usually very large and can effectively be taken as 
infinite. 



8 



massless, the lower endpoint corresponds to the case where energy- momentum conservation is 
satisfied by particles X and Y alone. On the other hand, the upper endpoint is obtained when the 
invariant mass of X and Y equals mx + my, which corresponds to the situation where X and Y 
are produced at rest in their overall center-of-mass frame. Thus, we have 



Eh 



El 



0, 



m , 



[mx -f 
2mA 



my 



(14) 
(15) 



For any fixed 7, the differential decay width in the energy of particle b in the laboratory frame 
is no longer a simple rectangle due to non-trivial h{Eh). For any specific laboratory frame energy 
Ef), contributions should be taken from all relevant values of E^ and weighted by h{Ef,)- This can 
be written as 



1 dT 

fdEb 



fixed 7 



dE,. 



h{Eh 



2EbV¥ 



(16) 



where 



E^ 



max 



mm 



El 



El 



Eh 



Eh 



7+ a/t^ 



Eh 



(17) 
(18) 



with Eh running from to E^^^ 



(7 + 



1 I . Again, since the visible particle is assumed 



massless, is zero and so the second equality in eq. (17) holds trivially. 

Finding an analytic expression for the location of the peak is difficult because of the model- 
dependence of h{Eh), and it follows that the precise location of the peak is also model-dependent. 
Nevertheless, we can still obtain a bound on the position of the peak for fixed 7. Suppose that 
we are interested in the functional value of the energy distribution at a certain value of Eh in the 
laboratory frame; according to the integral representation given above, the relevant contributions 
to this Eh come from a range of center of mass energies which go from E'^ to E"^', where these are 
defined by 



E'hil + - 1) = Eh; 
E'^{l-Vl^l) = Eh. 



(19) 
(20) 



Each energy contributes with weight described by h{Eh), as implied by eq. (16). 

and denote the corresponding energy in the laboratory frame 



Let us assume that E'^ 
as given by 



E> 



limit 



El 



X7-y^^). 



(21) 
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_ ^limit 

From these considerations, it follows that all rest-frame energies in the range from E', — 



to E'l^ = E^^^ contribute to a chosen energy in the laboratory frame, E']^™'*. On the other hand, any 

_ glimit _ _ 

laboratory frame energy greater than E} has contributions from E', > . to E'J = EV 



max. 



E. 



peak 
b 



the relevant range of the rest-frame energy values will shrink so that the peak cannot exceed El™'* : 

< ^""""{1 - - 1) < K'"'" for any fixed 7. (22) 

fixed 7 

In order to ensure that the first inequality holds even for 7 = 1, we assume in the last equation 
that h (E™^^) = 0, which is typically the case for a three-body decay. In order to obtain the shape 
of the energy distribution of particle h in the laboratory frame, all relevant values of 7 should be 
integrated over as with the two-body kinematics in the previous section. Hence, the laboratory 
frame distribution reads 

f3-hodyiEb) = -— = dEb _ d-f = . (23) 

I dEb Je< J^min{Et,E,) 2EfcV72-l 

Since the argument leading to eq. (22) holds for every 7, the superposition of contributions from 
all relevant boost factors does not alter this observation. Therefore, we can see that irrespective of 
^'(7) and h(Eb), the peak position of the energy distribution of particle b in the laboratory frame 
is always less than the maximum rest-frame energy: 

Ef ^'^ < Efe°^^^ . (24) 

To gain intuition on the magnitude of the typical difference between the peak of the energy 
distribution in the laboratory frame and the maximum rest frame energy, we show the ratio of 
the two as a function of 7 in Fig. 2. From the figure, it is clear that as the typical 7 increases 
beyond 7 = 1, i.e., as the system becomes more boosted, the location of the peak in the energy 
distribution becomes smaller. An appreciable shift of order 10% is achieved for a modest boost of 
order 7 — 1 ~ 10^^. 

It should be noted that all results here for both two-body and three-body decays are valid to 
leading order in perturbation theory. The presence of extra radiation in the decay will effectively 
add extra bodies to the relevant kinematics. Specifically, extra radiation can turn a two-body 
decay into a three-body one, which for our investigation would constitute a fake signal of two 
DM particles being produced in the decay of a heavy new physics particle. Therefore, we have to 
remark that in some cases, for instance, when the heavy new physics is typically produced with 
very small boost, the differences between the two scenarios of DM stabilization may be tiny and 
a study beyond leading order may be necessary. From Fig. 2 it seems, however, that the typical 
effect of the presence of two dark matter particles per decay of the heavy new particle is to easily 
induce an order one effect on the peak position. Therefore, we anticipate that such an effect would 
be much larger than the expected uncertainty from higher order corrections, which we estimate to 
be of order 10%. 
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7 



Figure 2: Relative separation of the peak of the laboratory energy distribution from the maximal 
energy in the center-of-mass frame of the three-body decay kinematics as per eq. (24). The hor- 
izontal red dashed line marks a 10% variation of the peak energy from the maximal value in the 
rest frame. 

Before closing this section, we emphasize that we shall use the right-hand sides of eqs. (7) 
and (24) as "reference" values to which the measurements of their respective left-hand side values 
(extracted from the energy distribution) are to be compared. In the next section, we show that 
such a reference value can, in fact, be extracted from an analysis of Mt2- 

2.2 The kinematic endpoint of the Mt2 distribution 

In this section, we review how the Mt2 variable is implemented for the two- and three-body decays 
of heavy particles produced at a collider. For our Mt2 analysis, we make further assumptions as 
follow: 

1) all massive decay products, i.e., particles X and Y in eqs. (4) and (5), are invisible; 

2) the mother particles A are produced in pairs; 

3) the entire decay process is symmetric in the sense that the mother particles are pair-produced 
and then decay to the same decay products, that is 

AA, A^Xh or A^hXY , (25) 

for the two-body decay and the three-body decay, respectively. 

The last assumption is especially relevant to make contact with the problem of distinguishing the 
Z2 and the Z3 dark matter interactions, as detailed in the introduction. 

2.2.1 Tw^o-body decay, one visible and one invisible 

The Mt2 variable generalizes the transverse mass to the cases where pair-produced mother particles 
each decay into visible particles along with missing particles (see Ref. [16] and references therein 
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for a detailed review). Specifically, it can be evaluated for each event by a minimization of the 
two transverse masses in each decay chain, under the constraint that the sum of all the transverse 
momenta of the visible and invisible particles vanishes. 

By construction, each of the transverse masses in both decay chains involve the mass of the 
invisible particle(s), and thus so does Mt2- Since a priori we are not aware of the invisible particles' 
masses, we are required to introduce a trial mass parameter into the definition of Mt2- We denote 
this trial mass by rh. The dependence of the definition of Mt2 on the trial mass makes it a function 
of fh. This function has been shown in Ref. [10] to have a kinematic endpoint 



M?^2l-body(^) = C^2-body + V^2-body + > (26) 

where the C parameter is given by 

2 2 

= 2mx • ^^^^ 
This C parameter can be deduced from eq. (26) by substituting the experimental value of the 
kinematic endpoint and the chosen trial DM mass. 

2.2.2 Three-body decay, one visible and tw^o invisibles 

As previously mentioned, for three-body decays we assume that the extra particle Y is also invisible. 
Therefore, as far as the detectable final state is concerned, the three-body decay looks like a two- 
body process. Since we are not a priori aware of the number of invisible particles involved in the 
decay process, a natural assumption is to hypothesize a single invisible particle per decay chain as 
in a two-body decay. In this context, we shall refer to this supposition as the "naive" Mt2 method 
(for three-body decay) [lU]. 

In each event, this three-body decay can be understood as a two-body decay process where the 
two invisible particles X and Y behave like a single invisible particle with an effective mass equal 
to the invariant mass of the system formed by particles X and Y . As is well-known, the invariant 
mass of the particles X and Y follows a distribution and ranges from mx + rny to ruA- Therefore, 
the overall kinematic endpoint in the corresponding Mt2 distribution arises when the invariant 
mass of the X-Y system is minimized [lU]. The theoretical expectation for -^'r2^3_body similar 
to that of the two-body decay: 



M?^2l-body(^) = ^^3-body + VC-l-body + > (28) 

where the C parameter is given by 

m\ - {mx + myf , . 

^^-^"'i^ ■ ^ ^ 

When comparing to the two-body case, two different features should be noted. First, given 

the same mother particle, visible state, and trial DM mass, the kinematic endpoint of the Mt2 
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distribution for the three-body process is expected to be smaller than that of the two-body process. 
This is because for the three-body decay, one more invisible particle, Y, is involved (see and compare 
eqs. (27) and (29), i.e., nix +'mY > mx)- Second, the fall-off of the distribution of the three-body 
process at the endpoint is faster than in the two-body process. This is because in the three-body 
case more kinematic constraints need to be satisfied to reach the kinematic endpoint [10, 11]. 

Before closing the Section, a further critical observation is in order. According to eqs. (26) 
and (28), we see that the observed values of M^f^ as a function of the various chosen trial DM 
masses (m) can be fitted with the same equation in both the two- and three-body cases: 

M^iX. = C+VC + m^, (30) 

where the parameter C can be extracted from the fit. This will be used in the following to extract 
the C parameter without making any assumption on the number of invisible products in the decay. 

The fact that the Mt2 endpoint can be described with the same parametrization in terms of 
a generic C parameter, as in eq. (30), is not surprising. In fact, for the two-body case in events 
near the endpoint each mother needs to have its decay products {b and X) emitted at the same 
rapidity (although the two mothers A can be at different rapidities) [iG]. Analogously for the 
three-body case, the two invisible decay products {X and Y) and the particle b produced at the 
same interaction vertex all need to share the same rapidity. In such a situation, the two invisible 
particles are kinematically equivalent to a single invisible particle, and so the decay can still be 
effectively reduced to a two-body decay. In this sense, M^2^ for the three-body case corresponds 
to the same kinematic configuration that gives the endpoint for the two-body case. However, it 
must be noted that the C parameter actually provides different information in the two cases. For 
two-body decays, the C parameter in eq. (27) is the same as the rest-frame energy of particle b in 
eq. (6), whereas for three-body decays, the C parameter in eq. (29) is the same as the maximum 
energy of particle b in the rest frame in eq. (15) 

{El for two-body decays 

(31) 
E'^^^ for three-body decays. 

This observation puts us in the position to extract the C parameter from the Mt2 distribution 
and compare it with the peak value in the energy distribution of the visible particle so as to test 
the nature of the decay. 

3 General Strategy to distinguish Z2 and Z3 

We now apply the above theoretical observation to the determination of the underlying DM sta- 
bilization symmetry. To pinpoint this stabilization symmetry, we study the energy distribution 

^Alternatively one can interpret the C parameter of the three-body decay as the analogy of the two-body case 
where the mass of the single DM particle is replaced by the mass of the effective single body made of the two DM, 
i.e. the sum of the mass of the two DM particles, as apparent from the comparison of eqs. (27) and (29). 
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of the particle b from the process defined in eq. (25). In particular, we exploit relation between 
this energy distribution and the distribution of the Mj'2 variable in the same process. As will be 
clear from the following analysis, the correlation between features of the distribution of these two 
observables will allow us to make a much firmer statement than merely utilizing one of them. 

In point of fact, the Mt2 distribution of the process eq. (25) could itself in principle be a good 
discriminator between Z2 and Z3 models. Indeed, as discussed in Section 2.2.2, the kinematic 
endpoint in the Mt2 distribution of the visible particles from a duplicate three-body decay, which 
is realized under Z3 symmetry, develops a longer tail than that of two-body decays, the latter being 
realized under Z2 symmetry. Therefore, a less sharp fall-off near the endpoint could be a sign of 
more than one invisible particle in the decay [10, 11]. However, shape analyses of the tail of the 
Mt2 distribution are rather delicate, especially in the presence of a background. Besides the issues 
raised by the backgrounds, there are also some inherent complications in using only the shape of 
the Mt2 distribution to determine the underlying stabilization symmetry. For example, the effects 
of spin correlation could change the shape of the Mt2 distribution, particularly the behavior near 
the upper endpoint of the distribution. In other words, a certain "choice" of spin correlation could 
alter the sharp edge of the Mt2 distribution in Z2 models, mimicking the typical distribution shape 
characteristic of Z3 models, and vice versa. 

Alternatively, one could try to use the energy distribution of the b particles in events from 
the process eq. (25). Recall that the distribution of the visible particle energy in their mother 
particle's rest frame is 5 function-like in Z2 models, whereas the distribution in Z3 models is non- 
trivial. Therefore, once the decay products are boosted to the laboratory frame from their mother 
particle's rest frame, the energy distribution for Z3 physics is expected to be relatively broader for 
a given mother particle. However, it is very hard to quantify the width of the resulting energy 
distributions in both Z2 and Z3 models because it is strongly model-dependent. In particular, the 
shape of the energy distribution in the laboratory frame is governed by the boost distributions of 
the mother particles, which are subject to uncertainties. Such uncertainties come from the fact 
that we are not a priori aware of the underlying dynamics governing the new physics involved in 
the process eq. (25), which affects, for instance, the production mechanism of the mother particles. 

In order to overcome the difficulties described above, we propose here a combined analysis of the 
two distributions. The goal is to obtain a more robust technique that is sensitive to the differences 
between the Z2 and the Z3 models but largely independent of the other details of the models. Also, 
we aim at formulating a method that is less demanding from an experimental standpoint and more 
stable against the inclusion of experimental errors. The analysis proceeds in two steps as explained 
in the following. 

From the data, one first produces the Mt2 distribution using a trial DM mass and extracts 
the kinematic endpoint M^J^^^^s • Then, by substituting the measured endpoint into the function 
given in eq. (30), one obtains the C parameter. As illustrated in eq. (31), the C parameter has 
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different physical implications depending on the stabilization symmetry of the DM. For the Z2 case, 
it is the energy of the visible particle in the rest frame of its mother particle, and by virtue 
of [14, 15], it is expected to be the value of the peak of the energy distribution in the laboratory 
frame. Alternatively, for a Z3 model the C parameter is an upper bound to the peak of the energy 
distribution in the laboratory frame. Therefore, the comparison between the extracted C parameter 
and the peak position in the b particle energy distribution enables us to determine whether the 
relevant physics is Z2 or Z3. This observation can be summarized as follows: 

for Z2 

for Z3. (32) 

Some remarks must be made about our proposal. First, the use of the distribution of M'T2 
is needed only to the extent that this is useful to extract the C parameter. In fact, in order 
to find the reference value needed for the comparison of eq. (32), any other observable that is 
sensitive to the relevant combination of masses could be used. Second, spin correlation effects 
do not change the location of the peak in the energy distribution of the b particle as long as the 
bottom partners are produced unpolarized, as discussed earlier. Additionally, although the overall 
shape near the endpoint of the Mt2 distribution could be affected by non-trivial spin correlation 
effects, the endpoint value is not. Furthermore, substantial errors in the determination of the Mt2 
endpoint can be tolerated. In fact, as shown in Fig. 2, the difference between the reference value 
and the typical peak of the energy distribution in a three-body decay is quite large. 

For the above reasons, we believe that compared with other methods which utilize only Mt2, 
the method presented here is more general and more robust in highlighting the different kinematic 
behavior inherent to the two different stabilization symmetries. 

In order to demonstrate the feasibility of the proposed analysis, we work out in detail an 
application of our method to the case of pair production of partners of the b quark that decay into 
a b quark and one or two invisible particles in the next section. 

4 Application to b quark partner decays 

In this Section, we study in detail the production of b quark partners, B' , and their subsequent 
decay into b quarks and one or two DM particles. As mentioned in the introduction, b quark 
partners occur in many well-motivated extensions to the SM. In the following, we apply the results 
of Sections 2 and 3 with the underlying goal of "counting" the number of DM particles in the above 
decay process. Although we employ DM and a b quark partner with specific spin for the purpose 
of illustrating our technique, we emphasize that our method can be applied for any appropriate 
choice of spins for the involved particles. In fact, the choice of spins does not alter our results so 
long as the mother particles are produced unpolarized. 
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Because the b quark partners are charged under QCD, the dominant production channel at 
hadron colhders would be via color gauge interactions, which guarantee that the b quark partners 
would be produced unpolarized and in pairs. Due to the fact that these particles are produced 
in pairs, the above results given for Mt2 are in force. Furthermore, the unpolarized production 
guarantees that the results of Section 2 can be applied to the energy distribution. 

In what follows, we consider the QCD pair production of heavy b quark partners at the LHC 
running at a center-of-mass energy y/s = 14 TeV, and we take as signal processes: 

pp B'B' bbxx for Z2 , (33) 

pp B'B' bbxxXX for , (34) 

where x is the DM particle. Once produced, we assume that each B' decays into a b quark and 
either one or two stable neutral weakly-interacting particles (see also Fig. 1). These processes will 
appear in the detector as jets from the two b quarks and missing transverse energy 

pp^bb + pT for both Z2 and Z3. (35) 

Note that our program is meant to be carried out only after the discovery of heavy b quark 
partner. In fact, our focus is not on discovery, but on determining what type of symmetry governs 
the associated decays of such a particle once the discovery is made, specifically in the bb+^ channel. 
In order to achieve this goal, a high integrated luminosity would be required to make a definitive 
determination of the underlying symmetry. Likewise, compared with the criteria necessary to claim 
the discovery of such a resonance, a different set of event selection conditions would be likely have 
to be used in order to make a definitive determination of the underlying stabilization symmetry. 

For our proof-of-concept example, we take m^/ = 800 GeV and = 100 GeV while noting 
that searches for scalar b quark partners such as Ref. [17] are in principle sensitive to our final state. 
Unfortunately, there is no available interpretation of this search in terms of a fermionic partner; a 
naive rescaling of the current limits on a scalar partner with mass of about 650 GeV shows that 
our choice of mass parameters might be on the verge of exclusion. However, we remark that our 
choice is only for the purpose of illustrating our technique, and can just as easily be applied to a 
heavier B' . 

There are several SM backgrounds that are also able to give the same detector signature as our 
signal. Since we require a double 6-tagging, the main backgrounds to our signal are the following 
three processes: i) Z + bb, where Z decays into two neutrinos, ii) + bb, where the W decay 
products are not detected, and iii) where again the two W's from the top decay go undetected ^. 
The first background is irreducible, while the latter two are reducible. 

To reduce these backgrounds to a level that allows clear extraction of the features of the b-jet 
energy and Mt2 distribution, we put constraints on the following observables: 

®By undetected we mean that the decay products do not pass our selection criteria or are legitimately undetected. 
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• PT,ji is the transverse momentum of the hardest jet in the event, 

• fjT = \ — YliPT.il is the missing transverse energy of the event and is computed summing over 
all reconstructed objects, 

• St = XT^A^ i^ transverse sphericity of the event. Due to the tendency of QCD to produce 
strongly directional events, the background processes typically have small sphericity, while 
decay products of a heavy B' are expected to be significantly more isotropic and hence will 
preferentially have a larger sphericity [18]. 

In general, the mismeasurement of the momenta of the observable objects used to compute ^ 
can produce an instrumental source of pr, as opposed to a "physical" source of ^ which originates 
from invisible particles carrying away momentum. The mismeasurement of ^ can grow as objects 
of larger px are found in an event, and it is therefore useful to compare the measured missing 
transverse energy with some measure of the global transverse momentum of the event. For this 
reason, we introduce the quantity ^ 

/ = ^/M^s where Meff = ^ + \pT \ + \pt ja I > 

which is expected to be small for events where the ^ comes from mismeasurements, but should 
be large for events where invisible particles carry away momentum. Furthermore, when the instru- 
mental originates mostly from the mismeasurement of a single object, the ^ is expected to 
point approximately in the direction of one of the visible momenta. Therefore, the events where 
the ^ is purely instrumental are expected to have a small 

A(/.(^,jets), 

which is the angle between the direction of the missing transverse momentum and any PTj- 
To select signal events and reject background events, we choose the following set of cuts: 



leptons with |r/;| < 2.5 and pri > 20 GeV for I = e, r, (36a) 

2 6-tagged jets with \7]b\ < 2.5 and prbi > 100 GeV, prba > ^0 GeV, (36b) 

> 300 GeV , (36c) 

St > 0.4 , (36d) 

/ > 0.3 , (36e) 

A0min(-^, bi) > 0.2 rad for all the selected 6-jets bi . (36f) 



Note that the our cuts are of the same sort used in experimental searches for new physics in final 
states with large pr, leptons and jets including 1 or more 6-jets (see, for instance, [19]). However, 

^Sometimes a slightly different quantity /' — ftf^^i \PT.i\ is used in the same context of our /. The two variables 
have the same meaning and give similar results. 
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Cut 


Zs (5 ^ 6x) 


Zs {B ^ 


Z + bb{Z ^ ui?) 


IVo mts 


159.75 


159.75 




Precuts 


139.89 


136.73 


2927 


pi^ > 100 GeV, pi^ > 40 GeV 


139.64 


133.76 


971.9 


^ > 300 GeV 


101.73 


69.01 


19.93 


/>0.3 


89.66 


65.21 


19.40 


A(/)min > 0.2 


88.95 


64.31 


18.81 


St > 0.4 


30.03 


16.07 


1.96 


2 6-tagged jets 


13.29 


7.18 


0.87 



Table 1: Cross-sections in fh of the signals and the dominant background Z + bh after the cuts 
of eqs. (36). The mass spectrum for the signals is m^' = 800 GeV and = 100 GeV. The line 
"No cuts" is for the inclusive cross-section of the signal. The line "precuts" gives the cross-section 
after the cuts > 60 GeV, > 30 GeV, rji, < 2.5, ARbb > 0.7 that are imposed solely to avoid a 
divergence in the leading order computation of the background. In the last line, the rate of tagging 
b quarks is assumed 66% [22]. 

notice that in our analysis, we privilege the strength of the signal over the statistical significance 
of the observation. As already mentioned, we imagine this investigation being carried out after the 
initial discovery of a B' has taken place. Hence, we favor enhancing the signal to better study the 
detailed properties of the interaction(s) of B' . For this reason, we cut more aggressively on ^ and 
St than in experimental searches and other phenomenological papers focusing on the discovery of 
B's (see, for example, [5]). 

We consider quarks separated by AR > 0.7 as jets. With this as our condition on jet recon- 
struction, the cuts of eq. (36) can be readily applied to the signals and to the Z + bb background; 
the resulting cross-sections are shown in Table 1. These cross-sections are computed from samples 
of events obtained using the Monte Carlo event generator MadGraphS vl.4.7 [20] and parton distri- 
bution functions CTEQ6L1 [21]. For the sake of completeness, we specify that in generating these 
event samples we assumed a fermionic B' and a weakly interacting scalar x- However, as already 
stressed, we anticipate that different choices of spin for these particles will not significantly affect 
our final result because the production via QCD gives rise to an effectively unpolarized sample of 
b quark partners. 

The estimate of the reducible backgrounds requires more work, as it is particularly important 
to accurately model the possible causes that make 

pp ^ tt ^ bb + X and pp — > + bb 

a background to our 26 -|- fJr signal. In fact, these processes have larger cross sections than Z + bb. 
However, they also typically give rise to extra leptons or extra jets with respect to our selection 
criteria in eq. (36). Therefore, in order for us to consider them as background events, it is necessary 
for the extra leptons or jets to fail our selection criteria. Accordingly, the relevant cross-section for 
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these processes is significantly reduced compared to the total. In fact, we find that tt and W bb 
are subdominant background sources compared to Z + bb. In what follows, we describe how we 
estimated the background rate from tt and W^bb. 

An accurate determination of the proportion of tt and W^bb background events that pass the 
cuts in eq. (36) depends on the finer details of the detector used to observe these events. However, 
the most important causes for the extra jets and leptons in the reducible backgrounds to fail our 
jet and lepton identification criteria can be understood at the matrix element level. We estimate 
the rate of the reducible backgrounds by requiring that at the matrix element level, a suitable 
number of final states from the tt and W + bb production fail the selections of eq. (36) for one of 
the following reasons: 

• the lepton or quark is too soft, i.e., pT,i < 20 GeV, pTj < 30 GeV 

• or the lepton or quark is not central, i.e. \'r]ij\ > 2.5 . 

Additionally, when any quark or lepton is too close to a 6 quark, we consider them as having 
been merged by the detector, and the resulting object is counted as a 6 quark (i.e., AR^i < 0.7, 
ARhj < 0.7), or if any light quark or lepton is too close to a light jet, they are likewise merged, and 
the resulting object is counted as a light quark (i.e., ARji < 0.7, ARjj < 0.7). In the latter case, 
the light "jet" resulting from a merger must then also satisfy the px and rj criteria given above for 
going undetected. 

Using our method to estimate the results on the backgrounds in Ref. [5] , the analysis of which 
was carried out with objects reconstructed at the detector level, we find that our estimates agree 
with Ref. [5] within a factor of two. Because we successfully captured the leading effect, we did not 
feel the necessity of pursuing detector simulations in our analysis. 

Estimating the reducible background after the selections in eq. (36), we find that tt and W + bb 
are subdominant compared to Z + bb. The suppression of the reducible backgrounds, and in 
particular, of tt, comes especially from the combination of the St and frp cuts. This is shown 
in Fig. 3, where we plot the ^ distributions of the three backgrounds under different St cuts: 
St > 0, St > 0.2, and the cut St > 0.4, which is used in our final analysis. Clearly, one can see 
that for a ^ as large as our requirement in eq. (36), the dominant background is Z + bb, and that 
in particular, the tt is significantly suppressed by simultaneously requiring a large ^ and moderate 
St cut (rightmost panel in the figure) . 

As the first step in our analysis, we compute the Mt2 distributions expected at the LHC for 
our two potential cases of new physics interactions, Z2 and Z3 . The distributions for the two 
cases are shown in Fig. 4. Since we found that with selections of eq. (36), the Z + bb process is 
the dominant background, as seen in the figure, we consider it the only background process. The 
two distributions have been computed assuming a trial mass m = GeV and have an endpoint at 
787.5 GeV and 750 GeV for the Z2 and the Z3 cases, respectively. Interpreting the distributions 
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Figure 3: ^ distributions for the three backgroun(is {Z + bb, + bb, and tt) with St cuts of 
increasing magnitude, St > 0.0, > 0.2, and > 0.4 from the left panel to the right panel. In each 
plot, the black solid, blue dotdashed, and red dashed curves represent Z + bab, + bb, and tt, 
respectively. 



under the naive assumption of one invisible particle per decay of the B' , we obtain from eq. (30) a 
C parameter that is 383.75 GeV and 375 GeV for Z2 and Z3 , respectively. These are the reference 
values that we need for the analysis of the energy distributions 

As the final step in our analysis, we need to compare the obtained reference values with the 
peaks of the energy distributions. These distributions are shown in Fig. 5. We clearly see that the 
location of the peak in the energy distribution the Z2 case coincides with the associated reference 
value, whereas for the Z3 case the peak is, as expected, at an energy less than the associated 
reference value. We remark that in the Z3 case, the peak of the energy distribution is significantly 
displaced with respect to the reference value. Therefore, we expect our test of the Z2 nature of the 
interactions of the B' to be quite robust under the inclusion of both experimental and theoretical 
uncertainties, such as the smearing of the peak due to the resolution on the jet energy, the errors 
on the extraction of the reference value obtained from the Mt2 analysis, and the shift of the peak 
that is expected due to radiative corrections to the leading order of the decay of the B' . 



5 Conclusions 

In this treatise, we studied the problem of the experimental determination of the general structure 
of the interactions of an extension to the SM that hosts collider-stable WIMPs. If these new 
particles are charged under a new symmetry and the SM particles are not, then the lightest such 
WIMP is stable and is concomitantly a candidate for the DM of the universe. In the context of such 

^We remark that as apparent from the figure, the signal rate is much larger than that of the background, and 
therefore the shape of the distribution expected at the LHC largely reflects the features of the signal. In this case, it 
seems particularly straightforward to extract the endpoint of the distribution. In other cases where the background 
is larger, the extraction of the endpoint may require a more elaborate procedure, especially for the Z3 case where the 
endpoint is much less sharp (see, for example, [23, 24, 10, 25]). 
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Figure 4: Mt2 distributions after the cuts of eq. (36). The chosen masses for the new particles 
are rriB' = 800 GeV and = 100 GeV. The left panel is for the Z2 signal while the right panel 
is Z3 (both in blue). In both cases, the background \s Z + hb (red). In both panels, the black line 
represents the sum of signal and background. The black vertical dashed lines denote the theoretical 
prediction for the endpoints. 
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Figure 5: Energy distributions of the b quarks after the cuts of eqs. (36). The chosen masses for 
the new particles are rriB' = 800 GeV and rriy- = 100 GeV. The left panel is for the Z2 signal, while 
the right panel is Z3 (both in blue). In both cases, the background is Z + hh (red). In both panels, 
the black line represents the sum of signal and background. The black vertical dashed lines denote 
the reference values extracted from the Mt2 distributions of Fig. 4 using eq. (30). 
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DM models, our work is thus relevant for the determination of the stabilization symmetry of this 
DM. In more detail, such models typically have heavier new particles that are charged under both 
the SM gauge group and the DM stabilization symmetry. Thus, these particles can be produced 
via the collision of SM particles, and will decay into DM plus SM particles. The number of DM 
particles in such a decay depends on the DM stabilization symmetry. Our goal was to devise a 
strategy to count this number of DM and thus probe the nature of this symmetry, based only on 
the visible part of the decays. 

To illustrate the technique, we studied models with fermionic b quark partners, i.e. colored 
fermions with electric charge —1/3 with sizable coupling to the b quark. In our example, we 
considered the case of b quark partners with mass at or below the TeV scale. The possibility of 
such is motivated by extensions to the SM that solve the Planck-weak hierarchy problem, since 
they contain top partners and, thus by SU{2)l symmetry, bottom partners. In the same model, 
it is also possible to have a WIMP DM. The b quark partners, as the typical states of the new 
physics sector, are charged under this stabilization symmetry and will then decay into a bottom 
quark, plus DM. Furthermore, thanks to their color gauge interactions, the b quark partners have 
a large production cross-section at hadronic colliders. Therefore the study of b quark partners is 
very well-suited to illustrate our technique. 

The literature on b quark partners thus far has only considered single DM in each decay chain, 
as would be the case in models where the DM is stabilized by a Z2 symmetry. However, in general, 
there can be more than one DM in this decay chain; for example, two DM are allowed in the 
case of a Z3 stabilization symmetry, albeit not in the case of a Z2 symmetry. So, the question we 
posed is whether we can distinguish the hypothesis of one vs. (say) two DM particles appearing in 
each of these decay chains. As mentioned above, in this way we can probe the nature of the DM 
stabilization symmetry. The question is non-trivial, because in either case the detectable particles 
produced are the same, and so is the signal of the b quark partners' production, i.e. bb + 

To distinguish between one and two DM in each b quark partner decay chain, the first result 
we used is that the measured Mt2 endpoints can be fitted by the formula eq. (30) irrespectively of 
how many DM particles are produced. The value of the free parameter obtained by fitting eq. (30) 
to the data is used in the next step of our analysis as follows. The second theoretical observation 
is that the peak of the distribution of the b quark energy in the laboratory frame is the same as 
the mother rest frame value for the two-body decay, but is smaller than the maximum value in 
the mother rest frame for the three-body decay. The crux is that the rest frame energy that is 
used as a reference value in this comparison is precisely the parameter obtained in the above Mt2 
analysis. Combining the above two facts, we showed that the peak of observed bottom-jet energy 
being smaller than (vs. same as) the reference value obtained from the Mt2 endpoint provides 
evidence for two (vs. one) DM particles in the decay of a 6 quark partner, and thus a Z3 symmetry 
can be distinguished from Z2. 
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We verified our theoretical observations in B' pair production and decay at the LHC. To assess 
the feasibility of the determination of the stabilization symmetry with our method, we simulated 
the signal and the dominant SM backgrounds. Using suitable cuts, we showed that the background 
in this case is due mostly to Z + bb. We studied in detail the case where the b quark partner has 
a mass niB' = 800 GeV and the invisible particles have a mass = 100 GeV. In this case, the 
background can be made small compared to the signal using the cuts of eq. (36). In Figures 4 and 5, 
we show the resulting Mt2 and b quark energy distributions relevant to our analysis. We observed 
that the peak in the b quark energy distribution for Z2 models is consistent with the reference 
value from the Mt2 endpoint, while that of Z3 models is apparently less than the corresponding 
reference value. The determinations of the peak of the energy distribution and of the reference 
value needed for our analysis are subject to uncertainties, e.g. those that propagate from the error 
in the determination of the Mt2 endpoint. However, the evidence for a Z3 stabilization symmetry 
comes from a difference between the peak of the energy distribution and the reference value. The 
theoretical prediction for this difference is large enough compared to the relevant uncertainties so 
that the proposed method seems to be quite robust, and should allow a clear discrimination of the 
stabilization symmetry of the DM. 

In future work we plan to extend the theory of Section 2 to deal with massive visible decay 
products. Thus, we shall be able devise a strategy to tell apart Z2 and Z3 stabilization symmetry 
in top quark partners decays into a top quark and invisible particles, which arise in the same 
scenario that we studied here. We also expect that our theoretical observation can be relevant in 
other applications, such as distinguishing two-body from three-body decays independently of the 
issue of DM. 
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